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Abstract 

■ 

| The 0(a s ) QCD corrections to the three-body decay width of the charged 

Higgs T(H + — >W + bb) are discussed in the MSSM model. Our calculations in- 
dicate that the standard QCD corrections to the three-body decay mode raise 
the width by about 12% and the supersymmetric QCD corrections (due to g, t, b 
| exchanges) can be comparable to or even larger than the standard QCD corre- 

ct ■ tions in some regions of the supersymmetric parameter space. This is mainly 

due to the effect of large left-right mixing of stop(i). It could significantly affect 



the phenomenology of the H + search. 



o 

Ph! I Introduction 

PS. 

Despite of the great success of the standard model (SM) in elementary particle 
physics, one important aspect, the Higgs mechanism, of the model has not yet been 
experimentally verified, and there is plenty of room for extensions. The SM possesses 
^ ■ the minimal Higgs structure of one doublet and only one neutral physical Higgs boson 

is left after SSB. One of the most interesting version of the extended SM is the 
minimal supersymmetric standard model(MSSM) which demands two Higgs doublets 
and predicts a charged Higgs pair(i/ ± ) Search for charged Higgs will give bounds 
on the parameter space of new physics models. Several groups have analyzed their 
experimental data and give lower bounds on the charged Higgs mass 0, 

The top quark decay provides a promising laboratory to look for Higgs boson 
because the top quark-Higgs Yukawa coupling is large. The H ± search strategy in top- 
quark decays has been based on the distinctive features of the channels t — > bH + — > 
b{r + u, cs), compared with the standard model decay t — > bW + — ► b(£v, ru, q'q), which 
assumes that the dominant decay channels of the charged Higgs below tb threshold 
are H + — ► t + v and H + — > csM. Recently, another important three-body decay 
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channel of the Higgs boson H + — > bbW + has been proposed by Ernest Ma, et ai.fj. 
They found that the above three-body channel is dominant for H + > lAOGeV and 
tan/? < 1 while the t + v mode dominates at large tan/5. The cs mode is relatively 
small at all tan/3 for H + > 140GeV . This can significantly affect the signature for 
top quark products. 

Accurate predictions of charged Higgs width including 0(a s ) corrections to the 
above processes are important to the phenomenology of the H + search. The one-loop 
effects (strong and electroweak) on the decay t — > H + b and t — > bW + have been 
discussed in detail in and respectively. The leading QCD corrections to the 
process H + — > cs is taken into account by substituting the quark mass parameters 
by the running masses at the H + mass scale Q which changes the charm quark mass 
from 1.5GeV to lGeV and considerably reduces the width of the H + —>■ cs process. 
Correspondingly, the leading QCD corrections to the above three-body decay channel 
should also be taken into account, which, so far as we know, is absent in the literatures. 
In this paper, we present the calculations of the QCD and SUSY-QCD corrections to 
the width of the three-body decay H + — > bbW + within the MSSM. We also calculate 
the one-loop corrections from charged Higgs particle which may be important for 
small tan/3. 

The paper is organized as follows. In section II, we present most of the formulae 
in our calculations and show some important analytical results. We also discuss the 
arrangement of the ultraviolet and infrared divergences. Some relevant formulae are 
given in Appendix A and B. In section III, numerical results and discussions are 
presented. 



II Analysis and Formulae 

We shall perform the calculations mainly in the MSSM. The related pieces of La- 
grangian are given in Appendix A. For simplicity, we neglect terms directly propor- 
tional to rrib but of course keep those terms which are singular in the vanishing mass 
limit (~ ln(m&)). In the final results after we include real gluon emission, such singu- 
larities are cancelled out. As a result of taking m& = in the MSSM, only left-handed 
bottom-squark(fek) enters our calculations while the left-right handed mixing of the 
top-squark (ii,,tji) is included in the calculations. 
First we define several quantities and symbols 

G = cot (3\V tb \\ (1) 
2m w 

A x = u{p 2 ),tP L v( Pl ), (2) 

A 2 = u{p 2 )i>zPLv{pi)p2 ■ e, (3) 

-43 = u(p 2 )$ 3 Plv(Pi)Pi ■ e, (4) 
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where cot ft = Vi/v 2 is the ratio of the vaccum expectation values of the two Higgs 
doublets, Vtb is the CKM matrix element, e is the polarization vector of W + boson 
and Vii Vii Vz are the four-vector momenta of b, b and W + respectively. 
The tree-level amplitude due to Fig. 1 is 



Mn = G 



jAi, 



1 2 1 ' 

Pi - mf 



(5) 



with p t =P2 + P3 and the corresponding width H is 



dTr 
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3g 4 mfcot 2 ft 



dsidsb 
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A m tv 



mf - st. 
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(s w - 2m 2 ) + 
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m h 
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(6) 



where sg, Sf, and sw are the 4-momentum squared transferred to the corresponding 
particles by H + Hj. 

The one-loop QCD and SUSY-QCD corrections to the process arise from the 
diagrams of Fig. 3 and Fig. 4, which involve virtual gluon and gluino exchanges. The 
self-energy corrections of external legs make no contribution to the width and we do 
not depict these diagrams. The amplitudes for those diagrams in Fig. 3 and Fig. 4 



arc 



M {9 ' 9) = G 



O'.s- 



1 



3tt (pi 



m 



M, 



(9,9) 



G 



O's 



Mi 9 ' §) = G 



3vr (p 2 t - mi) 



3vr (p 2 t - m 2 t ) 



17(9,9) A 

[F^~ 9) A l + F^~ a) A 2 ] 



(7) 
(8) 

(9) 
(10) 

where the indices g and g refer to the gluon and gluino corrections. — Fj 9 ' 9 ^ 

are functions of various kinematic invariants, which include the two-point, three-point 

A9,g) 17(9,9) 



M. 



(9,9) 



.a* 



G r±[ F (9>~9) Ai + F ^9> A2 + F K9,9) A 

3ir 



(9,9) 



^(9,9) 
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and four-point functions as defined in |L0 . We give the explicit forms of 

F y9,9,_pK 



in Appendix B. It can be seen from above that the amplitudes for the loop diagrams 
in Fig 4. have the same structures as their counterparts in Fig. 3. 

Now we turn to the mass and wave-function renormalization due to Fig. 2(a), (b). 
The self-energy parts can be written in the form 



^(9) = 



3tt 



2B (p 2 , m 2 q , A 2 ) - 2B 1 (p 2 , m 2 q , A 2 )) 



+ (45 (p 2 ,m 2 ,A 2 )-2)m g ], 
<9)rJ2\ p. , ^(g). 



-iE® = -itypf ( P 2 )P L + zf(p 2 )P R ) + £?V)], 



12) 



3 



^1\p 2 )Pl + ^f{p 2 )P R = ^[2(cos 2 9P L + sm 2 9P B )B 1 (p 2 ,ml,ml) 

+ 2(sm 2 6P L + cos 2 9P R )B 1 (p 2 , m% m%% (13) 

(y 

^V) = ^[m~ g sm29(B (p 2 ,mlml) - B (p 2 ,m P ,ml))}, ( 14 ) 

where we have introduced a gluon mass A to regularize the infrared divergences and 
9 is the mixing angle of left-right handed top-squark defined in Appendix A. 

We shall use the on-shell renormalization scheme throughout the paper [IT]. From 
(|TTD and (|T2|), we can easily derive the mass counterterms <5m 

5m (9) = ^m q [2B (ml m 2 q , 0) - 2£? 1 (m 2 ., m 2 ,, 0) - 1], (15) 

(y. 

+ ^ sin 29(B (m 2 q , m% m|) - B Q (m\, m 2 , m|))], (16) 

Tflq 

and the wave-function renormalization constants 

5Z i9) = ^.[l-2B (m 2 q ,ml\ 2 )-2B 1 (mlm 2 q ,\ 2 ) 

+4m 2 q (DB (m 2 q , m% A 2 ) - DB 1 {m\, m 2 q , A 2 ))], (17) 
8Zf = T,f + m 2 q {tf + tf)+2m q tf , (18) 
8Zf = Y,f + m 2 q {±f+±f) + 2m q ±f\ (19) 

with X = dX/dp 2 . In the above expressions we denote -Bx(X=0,l) as the two point 
functions and DBx as the derivative of Bx with respect to the momentum squared. 
We give the related two point, three point and four point functions in Appendix B. 

The counterterms for the coupling coefficients of the Higgs-top-bottom vertex and 
top-bottom-H^ + boson vertex are 

5li = -^+ 1 -(5Z^ + 5Z^), (20) 
^ = \(SZ^ + 6Z^), (21) 



Si!l = --^ + kszS + 5Zl% (22) 
Siw = ks2&+52&). (23) 



The amplitudes of the diagrams containing above counterterms will cancel out the 
ultraviolet divergences contained in the amplitudes for Fig. 3(b),(c) and Fig. 4(b),(c) 
completely. 
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To remove the ultraviolet divergence of the amplitude for Fig. 4(a), we introduce 
the following counterterm which corresponds to Fig. 1(c) 

iWtfgPL + 8Z%P R ) - mt tL + 2 0/jtR + 5mf}. (24) 

There is another way to get the same results. If we only substract 5m from the 
self-energy function S at the mass-shell, the remaining ultraviolet divergences will 
disappear automatically after we add all the amplitudes for loop diagrams (including 
the self-energy diagrams of the external legs) and the mass counterterms related with 
the H + tb vertex. 

The 0(a s ) contributions to the three-body decay width are given by the interfer- 
ence terms between higher order amplitudes and tree-level amplitude, i.e. 

5T= [ 2Re(Y, M (Qfl) M*), (25) 

J vhasesvace , „ „ 



e,c,cr 



where e represents the W-boson polarization and c, a represent the color and spin of 
b and b respectively. 

The above results still contain infrared divergences which are cancelled out by the 
infrared divergences in real gluon emission given by Fig. 5. In our calculations of 
the real gluon emission contributions, we closely follow the procedure discussed in 
where an energy cutoff AE is adopted to distinguish soft and hard gluon. The 
soft gluon contribution is calculated by soft gluon approximation and the hard gluon 



contribution is calculated by Monte Carlo methods [12|, |13[ . The contributions of soft 
gluon emission of Fig. 5(a),(b) are 

(26) 





a s j dT 


dsbdsb 


3tc 2 11 dsbdsi 


dr^ 


ot s j dT 


dsbdsb 


3n 2 2Z dsbdsi 



(27) 

respectively, while that of the interference term between Fig. 5(a) and Fig. 5(b) is 

dT {mter) = ga.^ dT 
dsbdsb 37r 2 dsbds^ 

where 

r d*k 2 PtPj 

v J 2u k { Pi k){ Pj ky iyj 

|k|<AS 

with ujk = Vk 2 + A 2 and AE the cutoff parameter. The explicit form of Iy is given 
in |T0]. The infrared divergence terms contained in ( p6|) (|27| ) are the same 

2a Sln 4AE 2 dV 
3n A 2 dsbdsb 
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The infrared divergence terms contained in (|28|) is 



4a s . 2 Pl -p 2 4AE 2 
— In — In 

6tt mt 



dTr 



A 2 dshdsh 



(31) 



The infrared divergences in 



) (|27|) will cancel out the the infrared divergences 
"" ! f, — f ; " ' will cancel out the infrared divergence coming from 
the contribution of the box diagram in Fig. 3. This can be seen from the analytical 



contained in SZ b y> and that in 



form of 5Z, 



(</) 



6Z, 



(<?) 



a, r . 1 , , A 2 , u 2 

^ _ _ + ln47r _ 7i; -2 In — -ln^-4, 

37T e m b m b 



(32) 



where D=4-2e is the space-time dimension, 7^ the Euler's constant and /x the 't Hooft 
mass parameter in the dimensional regularization scheme. The divergent part of the 
amplitude for the box diagram is contained in the folowing integral 



M 



IR 



on 



4-D 



d U k 



(k 2 - A 2 )((p 2 - k) - m 2 b )((p 1 + k) - m 2 b ){pj - mf) 



G 



a s -AmtPx ■P2u{pi)4Plv{p{, 



3tt 



Pt 2 - rn 2 



■C Q (ml,m 2 b , (p 1 +p 2 ) 2 ,m 2 b ,X 2 ,m 2 b ). 



(33) 



The analytical expression of Cq can be found in [[L4| . We can thus easily obtain the 
infrared part contained in the contribution due to the box diagram in Fig. 3. 



dY{box) 
dsf,ds b 



-4a, 



In 



mi 



In 



A 2 



dT 



3n 2pi ■ p 2 2pi ■ p 2 ds b ds b 



(34) 



It is evident that all the infrared divergences are cancelled out completely. 

At last, we consider the charged Higgs loop-corrections to the width because we 
are interested in the large cot (3 region where the coupling of charged Higgs to top 
quark and bottom quark is large. The corresponding diagrams are shown in Fig. 6. 
The substraction procedure is standard [jllj. We will discuss the numerical result for 
it in the next section. 



Ill Numerical Calculations and Discussions 

We now turn to the numerical evaluation of the corrected width. We have tested 
the results in a number of ways. We found that the results are reliable since they 
don't depend on the choice of the t'Hooft mass parameters \i in the dimensional- 
regularization scheme and the fictitious gluon mass A. We also choose different energy 
cutoffs in the calculation of soft gluon emission and find the results are independent 
of it to a satisfactory precision. The cancellation of In m& is examined too. 
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The whole analysis will depend on m#+, tan/3, A(— A t ),mg,mi L (— m~ hl ) and 
mi R . For simplicity, we have assumed mi L = m (fl in the calculations. We take 
\V tb \=l, m t =180GeV, m b =4.8GeV, m w =80.4GeV, sin 2 ^=0.23, a=l/128, a s = 
« s (150GeV)=0.113 and m^=91.2GeV. The SUSY parameters are constained to sat- 
isfy the lighter top squark >90GeV. 

The effects of the standard QCD and SUSY-QCD corrections of (8), (9), (10), 
(11) can be seen in Fig. 7, in which we take (m^, m g , A) = (200,400,-300,200)GeV. 
The results show the QCD corrctions raise the width by about 12%. The gluino 
corrections can raise the width by about 9% if we take the parameters listed above. 

In Fig. 8, we show the contour lines of ^f£^p in the A-/i plane for tan/3=l and 
(mH+,Trig, m^=mf L ) = (150,400,280)GeV. This correction has a strong dependence on 
fi. It changes sign near \x ~ OGeV. We can see for the chosen values of masses it can 
reach about 12% when A is near ±450GeV and \i is near ±450GeV. 

In Fig. 9, we show the dependence of <5I ^^" e ™°- > as a function of A and mi L for 
tan/3=l and (m H +, m g , /i) = (150,420,300)GeV. We can see the gluino correction can 
be about 15% when A is negative and m- th is about 100 GeV. 

In Fig. 10 we show a contour plot of ^f^^ in the //-m^ plane for tan/3 = 1 
and (m H +, m g , A) = (150,420,300)GeV. We find the corrections can reach 20% when 
H = -500 ~ -400 and m ih is about 100 GeV. 

In Fig. 1 1 we show a contour plot of ^f^™^ lYl the tan (3-m g plane for (m H + ,fJ,,rrit L , 
A) = (150, -300, 280, 300)GeV. Devided by tan/3=l.l, the left part of the graph is sim- 
ilar to the right part. The correction rises when tan (3 deviate from tan/3=l.l. Given 
a tan f3, the corrections increase up to about mg=650GeV and then decrease as m g 
inceases. 

The reason for the large contribution of 5T (gluino) is mainly because the vertex 
correction part of the gluino-exchange corrections is proportional to the H + tb coupling 
which can be enlarged greatly if the t-mixing parameters A and \i are large. 

In Fig. 12, we show the ratio of the three-body decay width including the standard 
QCD corrections to the width of the two-body decay (H + — > cs,t + u), in which we 
can clearly see the corrected three-body decay rises sharply with increasing m#+. It 
can reach more than 8 times the size of the two-body decay in the given area. 

In Fig. 13, we show the ratio of ST^\ the corrections due to gluino exchanges, to 
the width of the two-body decay (H + — > t + u), in which we take (m t ~ L , m g , //, A) = (200, 
400, — 300, 200)GeV. These curves indicate that the ratio tends to be smaller when 
tan/3 increases while it can be large when tan/5 is less than 1. 

Finally, we have calculated the dominant terms from the Higgs loop corrections in 
Fig. 6. We find the Higgs corrections are relatively small, which can lower the width 
by about 0.6% cot 2 j3 when m ff +=140GeV, and the ratio of the Higgs corrections 
to the tree-level contribution decreses when m#+ increases. Therefore, the Higgs 
corrections can be neglected if tan (5 is not too small. 
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In Fig. 14, we show the ratio of T + 5T^ + 5T^ H+ \the width including the 
standard QCD and Charged Higgs correcions) to the width of the two-body decay 
(H + — > t + v). If we measure the three channels H + — > t + u,H + — > cs and H + — > 
W + bb simultaneously and obtain tan (3 of the two-Higgs doublet model of type II from 
T C s/T T + u , we can check the validity of the two-Higgs doublet model of type II without 
supersymmetry by comparing the experimental data for T 3 /T T + U with theoretical 
results like those in Fig. 14. 

In summary, we have performed complete calculations of the 0(a s ) standard QCD 
and SUSY-QCD corrections to the width of H + — > bbW + . We have found that the 
QCD corrections raise the width by about 12% and the SUSY-QCD corrections can 
be comparable to or even larger than the standard QCD corrections and change signs 
as /i varies. This provides an effective way to distinguish the two-Higgs doublet model 
of type II from the MSSM. 
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Appendix A 



In this appendix we list some relevant pieces of the SUSY Lagrangian. 
The charged Higgs boson coupling to top and bottom quarks is given by 



gVtb 
V2mw 



H + [m t cot (3tP L b + m b tan f3bP R t] + h.c, (A.l) 



where Pl,r = ^(1 TTs) are the chiral projector operators, cot (3 = t>i/i>2 is the ratio of 
the vaccum expectation values of the two Higgs doublets and V tb is the CKM matrix 
element. 

The squark couplings to the charged Higgs, gluino and W + boson are given by 

C fbH = ^-[(m 2 w sin2/3-m t 2 cot/?-m2tan/3)4^ - 

+ m t (A* t cot (3 - fiZ^Zlj + m 6 (A b tan/3 - ^)Z\ 3 Z\^ h l-b)H' + h.c, 

(A.2) 

C m = g s V2t*T a ~g a [-Zl*P L + Z%P R ]t + h.c, (A.3) 
A* = g.VWT*r[-Z%PL + Z b 2 *P R ]b + h.c, (A.4) 
C Uw = ^Z^V^Tt^W- + h.c, (A.5) 

respectively, in which g s is the QCD coupling constant, g is the EW coupling constant 
and T a is the matrices of the SU(3) generators in the 3 representation. 

_ f / cos# — sin# \ , . „, 

Z = • a a A - 6 
y smfl cosy J 

and Z b are orthogonal matrices which diagonalize the mass matrices M~ and M~ of 
the squarks. 

M 2 = ( m l + m t + 0-35 (m| cos 2(3) -m t (A t + ^cot/3) \ 

* ^ -m t (A t + /i cot (3) m\ R + m 2 + 0.15 (m| cos 2/3) ) ' 1 J 

M 2 = / m? i +m 2 -0.42(m|cos2/5) -m b (A b + /itan/3) \ 

6 ^ -m b (A b + fi tan /3) m 2 ^ + m 2 , - 0.08(m| cos2/5) J ' 1 j 

when we neglect the mass of bottom-quark, the last mass matrix reduces to 



rn 



2 



0.42(m| cos 2(3) 



M i-{ bL 4 m? fl -0.08(m|cos2/?) J' (A ' 9) 
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Appendix B 



In below expressions, m b = m bi and Zij(i,j = 1,2) represent the elements of matrix 
Z l defined in Appendix A. 



F[ g) = -2m\ + 2m 3 t B (0,m 2 t ,0) + 2m t {m 2 + s b )B (s b , m 2 , 0), (B.l) 
F 2 {g) = -m t [l - 2B (m^,0, m 2 ) + 2(m 2 w - s b )Co(m 2 W) s b ,m 2 b ,m 2 b ,m 2 ,0) 

+ 2{m 2 w - 2s 5 )Ci(s s , m 2 w , m 2 b , 0, m 2 , m 2 b ) 

+ 2{m 2 w - s b )C 2 (s b , m 2 w , m 2 b , 0, m 2 , m 2 b ) 

+ ACoo(s b: m 2 V: m 2 b: 0,m 2 : ml)] : (B.2) 
F 3 (9) = 4m t [Ci(s5,m^,m^, 0,7^,77^) + C xx {s b ,m 2 w ,m 2 b ^,m 2 t ,m 2 b ) + 

Ci2{s b ,m 2 w ,m 2 b ^,m 2 t ,m 2 b )}, (B.3) 
F 4 (9) = -2m t [l - 2B (m 2 H+ ,0,m 2 ) 

- {m 2 w - s b - s w )C (s b ,m 2 H+ ,m 2 b ,0,m 2 ,m 2 b ) 

- {m 2 w + s b - s b - s H /)C 1 (s5,m^+,mfe,0,m^,m^) 

- {m 2 w - s b - s^)C2(s5,m^+,mfe,0,m t 2 ,m^)], (B.4) 
F^ 9) = 2m t C (m 2 v ,s w ,m 2 H+ ,m 2 ,m 2 b ,m 2 b ) 

- 2m t s w D (m 2 b , m 2 w , m 2 H+ ,m 2 b , s b , s w , 0, m 2 b , m 2 , m 2 b ) 

- 2m t s w D 1 (m 2 b , m 2 w , m 2 H+ ,m 2 b , s b , s w , 0, m 2 b , m 2 , m 2 b ) 

+ 2m t (s b -s b - s w )D 2 (m 2 b , m 2 W) m 2 H+ ,m 2 b , s b , s w , 0, m 2 b , m 2 , m 2 b ) 

- 2m t s w D 3 (m 2 b , m 2 W) m 2 H +,m 2 b , s b , s w , 0, m 2 b , m\, m 2 b ) 

- Am t D 00 (m 2 b , m 2 w , m 2 H+ ,m 2 b , s b , s w , 0, m 2 b , m 2 , m 2 b ), (B.5) 
F^ g) = -4m t D 2 {m 2 b ,m 2 v ,m 2 H+ ,m 2 b , s b , s w , , m 2 b , m 2 , m 2 b ) 

- Am t D 12 (m 2 b , m 2 w , m 2 H+ ,m 2 b , s b , s w , 0, m 2 b , m 2 , m\) 

- Am t D 22 (m 2 b , m 2 w , m 2 H+ ,m 2 b , s b , s w , 0, m 2 b , m 2 , m 2 b ), (B.6) 
F^ g) = AmtD^m^mw^H+^^s^sw^,™^™ 2 ,™^ 

+ Am t D 23 (m 2 b , m 2 W) m 2 H +,m 2 b , s b , s w , 0, m 2 b , m 2 , m 2 b ), (B.7) 
= m t B (0,mlml)(m 2 - m\^Z 2 u + Z 2 2l ) - B (s h m 2 ,ml)(m 2 m t Z 2 u 

+ mtSiZli - m t m 2 .Zl i - 2m g m 2 Z li Z 2i - 2mgS b Z X iZ 2i 

+ mlmtZ^ + mtSiZl - m t m\ Z 2i ), (B.8) 
= Am t [C m (mlm 2 w ,s b ,m 2 ,mlml)Z 2 u }, (B.9) 
F 3 {9) = AZ u \m t C l2 [m 2 b , m 2 W) s b , m], m\, m\)Z u 

+ m t C 22 (m 2 b , m 2 w , s b , m 2 g , m-, m\) 
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i rt I 2 2 2 2 2\ r? 

+ m g C (m w , si, m h ,m~ h ,m lo m g )Z 2i 
+ m~ g Ci{m 2 b , m 2 w , s b , rri~, m\, m~.)Z 2i 

+ C 2 (m 2 b ,m 2 w , s b ,m 2 ~,ml,m\)(m g Z 2i + m t Z u )}, (B.10) 
Fa ] = — — ^ m 2 H + , ml, m\, m\ , m\)Z u 
+ m g m t C (s b , m 2 H +,ml, m 2 g , m\,m 2 b )Z 2i ] 

[(mly sin 2/5 - m 2 cot + m t (A cot /? - //)Z 2i ], (B.ll) 

771(3) — ^ n/222 2 2 2 2 2\ 

^5 = ^D 00 {m b ,m w ,m H+ ,m b ,s- b ,s w ,m-m- b ,m i m- h ) 

m t cot p y 

Z u [(m^ sin 2/5 - m 2 cot + m t (A cot /3 - //)Z 2i ], (B.12) 

17(g) ~~ ^ rn / 2 2 2 2 222 2\ 

6 = ~ CQt p [ D 2{m b ,m w , m H+ ,m b , s b ,s w ,m g , m b ,m ii ,m b ) 

Zu(( m w sin 2/3 — m 2 cot (3)Z U + m t (A cot (3 — n)Z 2i ) 

+ D 12 (m 2 b , m 2 w , m 2 H+ ,m 2 b , s b , s w , m\, m\, m\,m\)Z u 

((mly sin 2(3 — m 2 cot /3)Z U + m t (A cot (3 — fj)Z 2i ) 

+ D 22 (m 2 b , m 2 w , m 2 H +,m 2 b , s b , s w ,m 2 g ,m 2 b , m\,m~)Z Xi 

((mly sin 2(3 - m 2 cot (3)Z U + m t (A cot (3 - fi)Z 2i )], (B.13) 
4 

m t cot /3 

^[(m 2 ^ sin 2/3 - m 2 t cot (3)Z U + m t (Acot/3 - /x)Z 2i ]. (B.14) 



771(9) ^n/222 2 2 2 2 2\ 

F 7 = 1 D 23 (m b ,m w , m H +,m b , s b , s w , m g , m^m^m-^ 



The relevant scalar fuctions are defined as follows 

B (pl, ml, m\) = (wr 2 )" W) 4 "* / d D q[(q 2 - m 2 )((q + Pl ) 2 - rn 2 )]" 1 , (B.f 5) 
Co(p 2 l ,Pi 2 ,p 2 2 ,ml,m 2 l ,m 2 2 ) 

= (tn 2 )-\2nf,) 4 - D Jd D q[(q 2 - m 2 )((q + Pl ) 2 - m 2 )((g + p 2 ) 2 - m 2 )]" 1 , 

(B.16) 

D (p 2 i,Pi2, P23,pI, P 2 2 , P13, ml,ml,m 2 2 ,ml) 
= ( l n 2 y\27ifi) 4 - D 

J d D q[(q 2 - ml)((q + Pl f - m 2 )((g + p 2 ) 2 - m 2 2 )((q + p 3 ) 2 - mj)}- 1 (B.l 7) 
nB / 2 2 2^ _ dB (p 2 ,m 2 ,m 2 ) 

VB {p 1 ,m ,m 1 ) - > (B.18J 

in which pjj = (p, — Pj) 2 ■ 

The definitions of the tensor-integrals and the relevant decomposions are given 
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below 

T n (p u ■ ■ ;Pn-i, mo, • • ; m*-i) = - [ d D q g ff ' - ^ , (B.20) 

with the denominator factors D = q 2 — m\, D{ = (q + pi) 2 — mf (i=l,- • -,N-1) 

B^Pi^Bu (B.21) 

C, = Pi.C, + p 2 ^C 2 = ^p^Q, (B.22) 

i=i 

CV = 9^Coo + P\^P\ V C X1 + p 2>J p 2v C 22 + {pi^P2 u + P2^Pi u )C 12 , (B.23) 

2 

= g^ c oo + J2 Pi^Piv C iv (B.24) 

3 

D^^Pi^Di, (B.25) 
i=i 

3 

= g^D 00 + ^ pi^p^Dij. (B.26) 
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Figure Captions 



FIG. 1 Tree-level diagram for the three-body decay of the charged Higgs H + — > W + bb. 

FIG. 2 Diagrams relevant for wave- function and mass renormalization in the cal- 
culations of the 0(a s ) QCD and SUSY-QCD corrections to the width of the 
three-body decay H + — > W + bb. 

FIG. 3 Diagrams relevant for the calculation of the 0(a s ) QCD corrections to the 
width of the three-body decay H + — > W + bb. 

FIG. 4 Diagrams relevant for the calculation of the 0(a s ) gluino corrections to the 
width of the three-body decay H + — > W + bb. 

FIG. 5 Diagrams for real-gluon emission H + — > W + bbg. 

FIG. 6 Diagrams relevant for the charged Higgs corrections to the width of the 
three-body decay H + — > W + bb. 

FIG. 7 Comparision of the tree-level width with those including the standard QCD 
or gluino corrections . The parameters taken in gluino corrections are (m^, m g , fi, 
A) = (200,400,-300,200). 

FIG. 8 Contour lines of 5T (gluino) /Y (tree) for (m H +, m g , m t J = (150, 400, 280)GeV 
and tan/3=l in the A-/i plane. 

FIG. 9 Contour lines of 5T (gluino) /V (tree) for (m H +, m g , A«)=(150, 420, 300)GeV 
and tan/3=l in the A-n%i h plane. 

FIG. 10 Contour lines of 5T (gluino) /T (tree) for (m H +, m 5 ,A) = (150, 420, 300)GeV 
and tan/5=l in the [i-m- tL plane. 

FIG. 11 Contour lines of 5T (gluino) /V (tree) for (m H +, //, m~ tL , A)=(150, -300, 280, 
300)GeV in the tan f3-m g plane. 

FIG. 12 The ratio of the three-body decay width including the standard QCD cor- 
rections to the width of the two-body decay (H + — > cs, t + u) for m^+=150GeV. 

FIG. 13 The ratio of the gluino corrections to the three-body decay width to the 
width of the two-body decay (H + — > t + u) for (m H +,mi L ,m g , /i, A) = (150, 200, 400, 
- 300,200)GeV. 

FIG. 14 The ratio of the the three-body decay width including the standard QCD 
and charged Higgs corrections to the width of the two-body decay (H + — > r + u) 
for m^+=150GeV. 
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